Correlation functions for diffusion-limited annihilation, A + A — > 
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The full hierarchy of multiple-point correlation functions for diffusion-limited annihilation, A + 
A —* 0, is obtained analytically and explicitly, following the method of intervals. In the long time 
asymptotic limit, the correlation functions of annihilation are identical to those of coalescence, 
A + A — > A, despite differences between the two models in other statistical measures, such as the 
interparticle distribution function. 
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The kinetics of nonequilibrium processes, in particular diffusion-limited reactions, have attracted much recent 
interest [1-8]. Because of the lack of a comprehensive approach for the study of such systems, models that yield 
to an exact analysis are of prime importance. In this respect, none have been studied more than diffusion-limited 
annihilation, A + A — > [9-28], and coalescence, A + A — > A [9,13-15,17,21,23,26,28-33]. Known exact results include 
the time dependence of the particle concentration and the two-point correlation function (for finding two particles at 
two different points, simultaneously). It has also been shown that the full hierarchy of n-point correlation functions 
for the two processes is identical [21,23,34,35], but explicit expressions for n > 3 are unavailable. 

Here we attack the problem of correlation functions for annihilation, using the method of parity intervals (or 
even/odd intervals) [20,26,27,36]. We recover the identity relation of the n-point correlation functions for annihilation 
and coalescence, and we derive explicit expressions, valid in the long time asymptotic limit, for all n. 

Consider the annihilation model, defined on the line — oo < x < oo. Particles A are represented by points which 
perform unbiased diffusion with a diffusion constant D. When two particles meet they annihilate instantly. Since the 
reaction step is infinitely fast, the system models the diffusion-limited annihilation process A + A — > 0. 

An exact treatment of the problem is possible through the method of parity intervals [20,26,27,36]. The key 
parameter is G(x,y;t) — the probability that the interval [x,y] contains an even number of particles at time t [37]. 
Particles near the edges of an interval may diffuse into or out of the interval, affecting the probability G. (On the 
other hand, reactions inside the interval cannot affect its parity.) With this observation in mind, one derives a rate 
equation for the probability G(x, y; t) [26]: 

^-G(x,y ) t)=D(^ + ^)G(x,y;t). (1) 

The annihilation reaction imposes the boundary condition 

lim G(x, y;t) = l, (2a) 

x"\y or y\x 



and G must also obey the conditions required from a probability density function. If the initial distribution of particles 

2 2 



is random, then G(x, y; 0) = \ + \ exp[— 2crj(y — x)], where cq is their initial density. In this case we have the additional 



boundary condition 



lim G{x, y-t) = \ , (2b) 

x — > — oo or y — >oo 



From G(x, y; t) one can derive the particle concentration: 



d_ 



p(x;t) = -—G(x,y',t)\ v=x . (3) 



Let p n {x\,X2, ■ ■ ■ ,x n ]t) be the n-point density correlation function, for finding particles at each of the locations 
X\,X2, ■ ■ ■ ,x n at time t. The particle concentration, p(x, t) = p\(x, t), represents merely the first term in the hierarchy 

{p n }, 71=1,2,.... 
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The correlation functions may be obtained from a generalization of the method of parity intervals, in the following 
way. Let H n (xi,yi, x 2 , y 2 , ■ ■ ■ ,x n ,y n ;t) be the joint probability that the interval [x±, y{\ contains an even number of 
particles, [£2,2/2] contains an odd number, etc., (odd intervals are denoted by an overbar), at time t. The intervals 
are non-overlapping, and ordered: x\ < y\ < ■ ■ ■ < x n < y n . Let F n (x\, yi, X2, J/2, • • • , x n , y n ; t) denote the probability 
that the total number of particles contained in 1J™ =1 yi], is even. Thus, 

Fi(xi,yi;t) = Hi(xi,yi;t) = G(x\,y\\t) , (4a) 
F2(xi,yi,x 2 ,y2]t) = H 2 (xi,yi,X2,y 2 ;t) + H 2 (xi,yi,x 2 ,y 2 ;t) , (4b) 

and, in general, F n is expressible as a sum of 2™ _1 H n functions, corresponding to the different combinations of 
interval parities that contribute to a total number of particles that is even. Then, in view of Eq. (3), the n-point 
correlation function is given by 



(-1)™ d n 

Pn (•£!> • • • ) x n\ t) = cy n ^\ ~~q^ q F n (x\, y\, . . . , X n , y n ] t)\y 1 =xi,...,y n =x n ■ (5) 



The H n satisfy a 2n-dimcnsional diffusion equation, analogous to Eq. (1), and for similar reasons. However, the 
boundary conditions of this equation are complicated by the following fact. For yi — * X{+i, a particle moving from 
the i-th interval to the (i + l)-th interval, or vice versa, would flip the parity of the two adjacent intervals. On the 
other hand, F n satisfies the same diffusion equation as H ni 

d d 2 d 2 d 2 d 2 

-F n (x 1 ,y 1 ,...,x n ,y n ;t) = D(^ + ^ + --- + — + —)F n , (6) 

but the boundary conditions are simpler: F n contains also the case where the parity of the intervals i and (i + 1) is 
flipped, so it is not affected by a particle hopping between the two intervals. If interval i is shrunk to zero, Eq. (2a) 
yields the boundary condition 

lim F n (xi,yi, . . .,x n ,y n ;t) = F„_i(xi, .,x n ,y n ;t) , (7a) 

Xilyi or yilxt 

where we use the notation that crossed out arguments (e.g. $i) have been removed. If the endpoints of two adjacent 
intervals are brought together, the intervals merge, resulting in the boundary condition 

lim F n (xi,y!, . . .,x n ,y n ;t) = F n _ 1 (x 1 ,y 1 , . . .,x n ,y n ;t) . (7b) 

yi1x i+1 or Xi+iiyi 

Finally, for a random initial distribution of particles, we have 

lim F n (x 1 ,y 1 ,...,x n ,y n ;t) = \ . (7c) 

x\ — > — oo or y n — >oo 

The F n are tied together in an hierarchical fashion through the boundary conditions (7a) and (7b): one must know 
F n -\ in order to compute F n . At the root of the hierarchy, F\ = G is obtained from Eqs. (1) and (2). 

The problem posed by Eqs. (1), (2), (6), (7) is similar to that of diffusion-limited coalescence, A + A — > A [32]. In 
that case one defines E n (xi,yi, . . . , x n , y n ; t) as the joint probability of finding the intervals \x%, yi], i — 1, 2, . . . ,n, 
empty at time t. E\(x, y; t) = E(x, y; t) satisfies the same equation as Eq. (1): 

^E(x, y;t) = D(^ + ^)E(x, y; t) , (8) 

with the boundary conditions 

lim E(x,y;t) = l, (9a) 

xly or y\.x 

lim E(x,y;t) = 0. (9b) 

x^ — oo or y^oc 

Note the difference between the boundary conditions (2b) and (9b). Likewise, E n satisfies the same equation as 
Eq. (6): 



d d 2 d 2 d 2 d 2 

^ t En(x 1 ,y 1 ,...,x n ,y n ;t) = D(^ + ^ + ... + — + —)E n , (10) 
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with boundary conditions analogous to (7a) and (7b), 

lim E n {x u y x , . . .,x n ,y n ;t)= E n _ 1 (x 1 ,y 1 , .,x n ,y n ;t) , (11a) 

Xi\yi or yilXi 

lim E n (x 1 ,y 1 , . . .,x n ,y n ;t)= E n ^ 1 (x 1 ,y 1 , ...,#», . . .,x n ,y n ;t) . (lib) 

yi\x i+1 or x i+1 lyi 



but 



lim E n (x 1 ,y 1 ,...,x n ,y n ;t) = , (11c) 

x\^ — oo or y n ^oo 



instead of (7c). The n-point correlation function for coalescence is 

pl oa \x 1 ,...,x n ;t) = (-1)" — ■ ■ ■ — E n (xi,yi, ...,x n , y n ; t)\ Vl = Xl ,..., Vn = Xn . (12) 

Eqs. (1), (2), (6), (7), and (8) - (11) imply that the solutions for F n and E n are simply related: 

F n (x\,y\, . . . ,x n ,y n ;t) = \ + \E n (x\,y\, . . .,x n ,y n ;t) , (13) 

provided that the same relation holds also for the initial conditions. Suppose that the initial distribution of particles 
is random, with initial concentration c annl for annihilation, and Cg oal for coalescence. Then E n (x\, . . . , y n ; 0) = 

exp[-c|f al (yi - xi H h y n - x n )], while F„(a;i, . . . , y n ; 0) = \ + \ exp[-2c anni (yi — x\-\ V y n - x n )]. Thus, the 

relation (13) is satisfied if Cq 11111 = ic§ oal . Moreover, Eqs. (5) and (12) imply that in this case 

pT ni (x 1 ,...,x n] t) = ^r P c n oa \x 1 ,...,x n] t) , (14) 

for all n. In other words, the n-point correlation functions for annihilation and coalescence are identical, as already 
found by others [21,23,34,35]. 

We now produce explicit expressions for the n-point correlation functions in the long time asymptotic limit. Recall 
first the solution for E n . For n = 2 the solution is [38], 

E 2 {xi,yi,x 2 ,y 2 ;t) = 

E(x 1 ,y 1 ;t)E(x 2 ,y 2 ;t) - E(x 1 ,x 2 ;t)E(y 1 ,y 2 ;t) + E(x 1 ,y 2 ;t)E(y 1 ,x 2 ;t) , (15) 
where E(x,y;t) is the solution of Eqs. (8), (9). Generally, for n > 2 [32], 

(2ti-1)!! 

t)= ,pi z i,p\ t) ■ ■ ■ E(z 2n -l tP , Z2n,p! t) , (16) 

P =i 

where Zi iP , z 2 , p , . . . , z 2n ^ p is an ordered permutation, p, of the variables x\, yi, . . . , x n , y n , such that 

Zl,p < Z 2 . p , Z^.p < Z4,p, ■ • ■ , z 2n-l,p < z 2n,p, and < Z3 iP < Z5 jP • • • < Z2n-l,p • (17) 

There are exactly (2n — 1)!! = 1-3 (2n — 1) such permutations. a p is +1 for even permutations (permutations 

that require an even number of exchanges between pairs of variables), or —1 for odd permutations. Alternatively, the 
E n may be obtained through the recursion relation: 

n 

E n (xi , 7/1, ... , x n: y n ] t) = + E{x l ,y j ;t)E n - l {]k l ,y l , . . . , Xj , fa . . . , x n ,y n ; t) 

3 = 1 

n 

-'Y / E(x 1 ,x j ;t)E n _ 1 (}6 1 ,y 1 , f6j,yj, . . .,x n ,y n ;t) , (18) 

3=2 

then p n may be computed through the relation (12), or (5) and (13). 
Consider the long-time asymptotic limit, where 

E(x, y] t)=etfc(-^L). (19) 
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Then, the long-time asymptotic n-point correlation function is: 

(2ra-l)N 

(Xi 

j 2/1) • • • ) ^n; Urn 

t) , (20a) 



t—>oc 

p=i 



where 



_ ( x -f) — i l/V87r_Dt annihilation, 
\ l/V^nDt coalescence, 



(20b) 



C(z u z 2 ;t) 



-1 (^1,22) = {xk,Vk) , 

erfc(£/ fc ) (zi, z 2 ) = (x k ,xi) , 

(«i,a&) = (a*,w), (20c) 



-e 



(^i, z 2 ) = (2/fe,a;0 , 
> --A&fce"^ (21, z 2 ) = {Vk,yi) , 



and we used the notation = (xi — x^)/ VSDt. For example, for n = 2,3, we get the long-time asymptotic 
expressions: 



- - ■ , ..s,,. 'erfcfei) , (21a) 



p 3 



+Vn&i(e ?21 - e * 3 i)erfc(5: 



21; 



+^62(6 ? ' 2 -e & ? 3i)erfc(^3 2 ) 

+\/7i : C3i(e"^ 1 - e _ &-&)erfc(£ 3 i) . (21b) 

In summary, we have confirmed the fact that the infinite hierarchies of n-point correlation functions for coalescence 
and annihilation are identical, using the method of parity intervals. The simplicity of our approach allowed us to 
obtain explicit expressions for the long-time asymptotic limit, given in Eqs. (20). We note that our results arc not 
restricted to long times. Indeed, for the case of a random distribution of particles, such that the initial concentration 
for annihilation is half that of coalescence, the identity holds at all times. In this case, explicit expressions for the 
correlation functions (valid at all times) can be obtained by using the full solution of Eqs. (8), (9), 

E(x,y;t) = erfc(^) - \ e ^ {e ^[l - erf( ' ~ ^_fP** )] - e-<«<->[l + crf( ^^^ )]} , (22) 

instead of the asymptotic expression of Eq. (19). 

Remarkably, the particle distributions in coalescence and annihilation differ, despite the correspondence of the 
n-point correlation functions. The probability density function, p(x), for the distance x between two neighboring 
particles illustrates this difference. For large x, p{x) ~ e~ x for coalescence, while p{x) ~ e~ x for annihilation. 
Evidently, the complete hierarchy of n-point correlation functions is not sufficient to determine an infinite-particle 
system uniquely, and p(x) cannot be computed from a knowledge of the p n . A study of p(x) and {p n } in finite systems 
might illuminate this curious phenomenon. 
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